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I Abstract. We discuss experimental features of identified particle production from nucleus- 

nucleus collisions. These features reflect hadronization from a deconfined partonic matter whose 
particle formation scheme is distinctly different from fragmentation phenomenology in elementary 
collisions. Multi-parton dynamics, such as quark coalescences or recombinations, appear to be es- 
sential to explain the experimental measurements at the intermediate transverse momentum of 2-5 
GeV/c. Constituent quarks seem to be the dominant degrees of freedom at hadronization. Heavy 
quark production should help quantify deconfined matter properties. 
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The advent of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab- 
oratory (BNL) has turned the next page in search for, and study of, the Quark Gluon 
Plasma (QGP). Most recent dAu baseline reaction measurements by all four RHIC ex- 
periments [1, 2, 3, 4] confirmed that a dense strongly interacting medium has been cre- 
ated in central AuAu collisions at RHIC. The QCD nature of the dense matter created 
at RHIC and whether the current experimental evidence proves the discovery of the 
QGP have been since under debate within the heavy ion physics community [5, 6, 7, 8]. 
We will address here primarily the physics 'soft' hadron experimental results and the 
related evidence for deconfinement. The major topics we address are: 1) features of 
hadronization and other evidence for a color deconfined bulk partonic matter; 2) the 
QCD properties of the matter at hadronization; and 3)charm production and related fu- 
ture experimental measurements capable to further quantify the properties of deconfined 
phase. 

Among topics we discuss in depth are: mechanisms of hadronization, experimental 
status of strange particle production enhancement, azimuthal parton flow anisotropy, 
bulk properties of dense hadron matter at parton fireball breakup, charm experimental 
status at RHIC. 



This report constitutes the combined contribution of the authors to the "VI Quark Confinement and 
the Hadron Spectrum" conference, held 21-25 September 2004 in Sardinia, Italy covering their oral 
presentations and the 'Deconfinement' discussion session; proceedings to be published by the American 
Institute of Physics. 
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2. FEATURES OF HADRONIZATION OBSERVED AT RHIC 



2.1. What we know about hadronization 

We observe, in the final state, a multitude of hadrons, irrespective of what happened 
and which reaction system is observed. The paradigm emerged that QCD color charges 
are confined and hadrons exist in color singlet state. Individual reactions occur between 
leptons and/or quarks. However, the conversion, i.e., hadronization of quarks and gluons 
(partons) has not been understood based on first principles. There are several widely 
studied models: 

1. Hadron formation in elementary and qq (that is nucleon-nucleon) collisions 
has been described in the pQCD domain in terms of several components. The 
particle production process is factorized into parton distribution functions, parton 
interaction processes and in final step, fragmentation functions for hadron forma- 
tion. 

The fragmentation function is assumed to be universal and can be obtained phe- 
nomenologically from e + e~ collision measurments. A typical Feynmen-Field [9] 
fragmentation process involves a leading parton of momentum p, which frag- 
ments into a hadron of momentum p^ whose properties are mostly determined 
by the leading parton. The fragmentation function is a function of variable 
z = Phi Pi z G (0,1) . Note that baryon production is found to be significantly 
suppressed compared to the production of light mesons: the baryon to pion ratio 
increases with z, but never exceeds 20% [10], see also the pp-STAR results in 
figure 2 below. 

2. In the soft (low pj_) particle production region, the pQCD framework and factor- 
ization break down, particles are believed not to be from fragmentation of partons. 
In the elementary interactions physics, string fragmentation models were inspired 
by the QCD description of the quark and anti-quark interaction. The Lund string 
model is one of the popular hadron formation models which has been successfully 
implemented in Monte Carlo description of e + e ~, pp, and nuclear collisions [11]. 
Note that, in the string fragmentation models, the baryon production is also sup- 
pressed because baryon formation requires the clustering of three quarks [12]. 

3. Soft multi particle production in particular in the AA (nuclear) collision domain is 
described within the (Fermi) Statistical Hadronization (SH) model. SH is a model 
of particle production in which the birth process of each particle fully saturates 
(maximizes) the quantum mechanical probability amplitude, and thus, all hadron 
yields are determined by the appropriate integrals of the accessible phase space. 
The statistical hadronization model introduced in 1950 [13, 14] has matured today 
to a full fledged tool for soft strongly interacting particle production, capable to 
describe in detail hadron abundances once the mass spectrum of hadron resonances 
is included [15]. The key SH parameters within the grand canonical formulation of 
particle phase space are the temperature T and, at a finite baryon density present 
in a AA reaction system, the baryochemical potential /is- It is generally accepted 
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that as the energy of the colliding nuclei varies, a wide domain of T and /1b 
is explored, see figure2 in [16]. Note that the ratio of baryons to mesons is, in 

'"baryon 

chemical equilibrium, also suppressed, since e t <1 

4. The reaction picture of a soup of quarks in a dense expanding fireball inevitably 
triggers development of recombination and coalescence models [17, 18]. Re- 
cent experimental RHIC results triggered further developments on quark coales- 
cence [19, 20] and recombination [21, 22] which all have the essence of multi- 
parton dynamics for hadron formation, despite significant differences in details. 

2.2. (Anti)Baryon yield 

Strangeness plays a particularly important role as a characteristic QGP observable. 
The enhanced production of (strange) (anti)baryons has been an expected feature in re- 
combinant hadronization of QGR In nuclear collisions, the strange hadron yield derives 
from two independent reaction steps following each other in time: 

1. the establishment of a ready supply of strange quark pairs s, s which occurs pre- 
dominantly in the initial hot phase of bulk partonic matter by the thermal pQCD 
gluon fusion processes gg — > ss [23], in a manner independent of the production of 
final state hadrons; 

2. the high initial s, s yield survives the process of fireball expansion evolution [17], 
and contributes in hadronization of pre-formed s, s quarks to an unusually high 
multistrange (anti)baryon abundance [24]. 

This is illustrated in figure 1. 




FIGURE 1. Illustration of the two step mechanism of strange hadron formation from QGP: inserts show 
gluon fusion into strangeness, followed by QGP recombinant hadronization. 
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Experiments indeed show very significant new features of hadron formation in AA 
interactions compared to elementary e + e~ and qq (that is pp) collisions. 

• Baryon production from AA collisions, especially for multi-strange hyperons, has 
been measured to be much larger than theoretical (superposition, cascade) model 
calculations. The hyperon production per number of participant pairs from AA 
collisions at the SPS is significantly enhanced in comparison with the value from 
pp collisions. The enhancement factor increases from A to Q. hyperons and with 
the collision centrality [25, 26]. 

• The increase in the baryon production from AA collisions has become much more 
prominent at RHIC energies. Figure 2 shows the ratios of p/n and A/Ks from 
central AuAu collisions at v /*aw = 130, and 200 GeV measured by PHENIX [27] 
and STAR [28, 29, 30]. 
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FIGURE 2. Ratios of A to Kg from AuAu and pp collisions (STAR) and ~p to n from AuAu collisions 
(PHENIX) as a function of transverse momentum (p±). In addition to resonance contributions in all 
hadrons, the A includes contributions from E° decays. 



The apparent difference, in figure 2, between the STAR and PHENIX ratios can be 
understood since the STAR A data include the electromagnetic decay contribution from 
E° and both comprise different 'towers' of hadron resonance decays. Note also that 
the PHENIX data are corrected for post-reaction weak decay contributions Y,Y — > AW. 
Similarly, the STAR A data have been corrected for feed-down contributions from multi- 
strange hyperon decays. Some early data, for example [31], are not shown because these 
data were not corrected for weak decay contributions. 

The large baryon to meson ratio, seen in figure 2, cannot be accommodated by the 
traditional, e.g., fragmentation scheme. The large ratio at the intermediate p± region 
provide clear evidence that particle formation dynamics in AA collisions at RHIC are 
distinctly different from the traditional hadron formation mechanism via fragmentation 
processes developed for the elementary e + e~ and nucleon-nucleon collisions. 
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The recombination models [21, 22] provided a satisfactory description of the particle 
yields, in particular, the large production of baryons in the intermediate p± region. The 
formation of a dense partonic system provides a parton density dependent increase in 
baryon yield as a function of collision centrality through the coalescence mechanism. 

Comparing AuAu with pp reactions at = 200 GeV, see figure3, we see another 

large baryon production enhancement for strange (anti)baryons A, A, reported 
by STAR [32]. This enhancement increases with centrality [33, 34] and with greater 
strangeness content as found in strangeness recombination model [17]. These feature 
were recognized early on as a characteristic signature of QGP [24] . 



CD 
> 

P10 



T3 



STAR 

Vs NN =200 GeV 



• A 



T 
■ 

1 



I if 




Cl 
Cl 



CD 
> 

"co 

110 



T3 



STAR 

Vs NN =200 GeV 



• A 
■ =+ 



10 



part 



; 1 ■ 



(b) 



10 



N 



part 



FIGURE 3. The STAR yields per participant for production of A and E~ on left (a) and A and S + on 
right (b) in AuAu collisions at -^/snn = 200 GeV. Error bars are statistical. Ranges for pp reference data 
indicate the systematic uncertainty. 

This pattern of baryon enhancement has been observed by the WA97 [25] and NA57 
[26] experiments at lower reaction energy available at the CERN-SPS. There is a gradual 
increase in the strange antibaryon yield with reaction energy [26] . The large difference 
between baryon and antibaryon yields at SPS energy range are due to the presence of 
a significant baryon density at lower reaction energy. For purpose of comparison with 
RHIC results, this baryon density effect can be removed by considering the geometric 
mean of the baryon and antibaryon yield. 

In our opinion, the results presented above for the systematics of strange hadron en- 
hancement demonstrate that same novel mechanism operates in the entire collision en- 
ergy interval spanned by these data. Only a deconfined quark-gluon plasma is a copi- 
ous source of these often rarely produced hadrons. When the plasma fireball breaks up 
into final hadrons, the high abundance of strange quarks and antiquarks manifests itself 
yielding high abundances of multi strange hadrons. In conventional reaction schemes, 
the production of particles containing two or more strange quarks is suppressed by the 
rarity of the required reaction processes. 
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2.3. Nuclear modification factor 



The nuclear modification factor is defined in a collision number scaled comparison of 
peripheral with central head-on collisions: 



Rcp 



[yield/N bin 



central 



[yield /N bin ]P eri P heral ' 
The nuclear modification factor has also been defined by, 

[yield]** 



Raa 



N bin x\yield]PP' 



(1) 



(2) 



where Nbm is the number of binary nucleon-nucleon collisions. The [yield]** and 
[yield] pp are particle yields (d 2 n/dp±dy) from AA and pp collisions, respectively. 

A Raa or Rcp of unity implies that particle production from AA collision is equivalent 
to a superposition of independent nucleon-nucleon collisions. Hard scattering processes 
within the pQCD framework are believed to follow approximately the binary scaling in 
the kinematic region where the nuclear shadowing of parton distribution function and 
the Cronin effect are not significant. 

Measurements of charged hadrons and neutral pions have revealed a strong suppres- 
sion at high p± region in central AuAu collisions [35, 36, 37, 38]. Recent dAu measure- 
ments [1, 2, 3, 4] have demonstrated that the large suppression of high p± particles in 
central AuAu collisions is mainly due to energy loss, presumably of partons traversing 
the dense matter created in these collisions. 

Au Au Vs NN = 200 GeV 
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FIGURE 4. Rcp of K ± , K$, K*, 0, A, E and £2 in comparison with charged hadron in dashed line. 
Distinct meson and baryon groups are observed. 

Figure 4 shows the Rcp of Kg, K*, 0, A, S and Q. as a function of p± from AuAu 
collisions at ^/snn = 200 GeV measured by the STAR collaboration, where the ratio is 

derived from the most central 5% to the peripheral 40 60% collisions. The dashed 

line is the Rcp of charged hadrons for reference. In the low p± region, soft particle 
production is dominated by the number of participant scaling. 
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The particle type dependence can be described using hydrodynamic flow which pre- 
dicts a mass dependence for the low p± spectra. In the intermediate p± region of 2 to 
5.5 GeV/c the p± dependence of Rep falls into two groups, one for mesons and one for 
baryons. Despite the large mass differences between ^5 and K*/(j), and between A and 
S little difference among the mesons and among the baryons has been observed within 
statistical errors. Particle dependence in the nuclear modification factor disappears only 
above a p± of 6 GeV/c, consistent with the expectation from conventional fragmenta- 
tion processes. The unique meson and baryon dependence in the intermediate p± region 
indicates the onset of a production dynamics very different from both fragmentation at 
high p± and hydrodynamic behavior at low p±. 



2.4. Azimuthal anisotropy 



The azimuthal angular particle distribution can be described by a Fourier expansion 

d 2 n 



p±dp±d(j) 



(l+2£v n cosn(0 -¥«)), 



(3) 



where is the reaction plane angle and V2 has been called elliptic flow [39]. In a 
non-central AA collision the overlapping participants form an almond shaped particle 
emission source, see figure 5. The reaction plane is defined by the vectors x (impact pa- 
rameter direction) and z (beam direction). The pressure gradient is greater and particles 
would experience larger expansion force along the short axis direction (in plane) than 
along the long axis (out plane), see figure 5 left, resulting in the final state in an ellipsoid 
in transverse momentum space. Theoretical, (typically hydrodynamic) model calcula- 
tions find that this p± dependent component in V2 is generated mostly at the early stage 
of the nuclear collision. 
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FIGURE 5. Schematic diagram to show two dynamical origins of angular anisotropy, one from hydro- 
dynamical expansion and the other from surface geometrical phase space. 

There is another dynamical mechanisms responsible for V2 due to the spatial geomet- 
rical component in the phase space of emitting source. Particle production from a dense 
matter can be squeezed more in the reaction plane than that out of the reaction plane. 
More generally, a freeze-out geometry will impose itself on particle abundance elliptic- 
ity. However, any p± dependence would be a 2nd order effect associated e.g. with the 
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optical depth of the freeze-out system. Figure 5 shows a schematic diagram for these two 
mechanisms of generating azimuthal angular anisotropy V2 in non-central collisions. 

Figure 6 shows the V2 as a function of p± for it, K, p, A and S from PHENIX [40] 
and STAR [41] measurements. The angular anisotropy V2 reveals three salient features: 

1 . particles exhibit hydrodynamic behavior in the low p± region — a common expan- 
sion velocity may be established from the pressure of the system and the heavier 
the particle the larger the momentum from the hydrodynamic motion leading to a 
decreasing ordering of V2 from 7t, to K and p for a given p^; 

2. V2 values do not depend strongly on p± at the intermediate p± region in contrast to 
the strong p± dependence in the yield of particles; 

3. the saturated V2 values for baryons are higher than those for mesons and there is a 
distinct grouping among mesons and baryons. 
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FIGURE 6. Azimuthal angular anisotropy V2 as a function of p± for identified particles. The hydrody- 
namic calculation results are by P. Huovinen et al. [42]. 

The absence of a strong p^ dependence at the intermediate p± region is an intriguing 
phenomenon for the angular anisotropy of particle emission. Parton energy loss in 
the dense medium created in AA collisions was proposed as a possible mechanism 
for generating an angular anisotropy v%. High energy partons are quenched inside the 
dense medium and lead to an effective particle emission just from a shell area of the 
participating volume [43]. Jet quenching scenario [44] (parton energy loss scenario) 
cannot explain the particle dependence in both the nuclear modification factor and the 
angular anisotropy V2 at the intermediate p± region. Within the parton energy loss 
scenario the larger V2 of baryons implies a higher energy loss than that of mesons; the 
larger nuclear modification factor of baryons, however, is only consistent with a smaller 
energy loss than that of mesons. This scenario may be important for the considerable 
V2 magnitude for charged hadrons at a p± greater than 6 GeV/c or so measured by 
STAR [45]. 
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Indeed, the magnitude of the measured V2 is significantly larger than what can be 
accommodated based on particle emission from a geometrical ellipsoid source within 
an energy loss scenario [46]. Using a more realistic Wood-Saxon description of the 
colliding nuclei for the ellipsoid source, the predicted theoretical V2 is much smaller than 
that from a hard-sphere model of the colliding nuclei, leading to a greater discrepancy 
between the measurement and the theoretical expectation. 

We conclude that the magnitude and the particle dependence of V2 at the intermediate 
p± region cannot have a dynamical origin either from hydrodynamic flow or from 
parton energy loss alone. Another plausible procedure would be to relate the V2 in the 
intermediate p± region to the geometrical shape of the emitting particle source, see 
figure 5. A surface emission pattern from the almond shaped participant volume should 
naturally lead to a saturation of vi- A surface emission scenario is possible if particles 
are produced by surface related dynamical instabilities and the hadronization duration is 
relatively short [47]. 



2.5. Quark Flow Anisotropy 

An empirical quark number n scaling of V2 has been noted [48]. n is the number of 
valance quarks and antiquarks in a hadron. Figure 7 presents the V2/n as a function of 
p±/n for K, p, A and S from AuAu200 collisions, where the line is a polynomial fit 
to the data points. The bottom panel shows the ratio of data points to the fit. At the 
intermediate p± region ( 0.6 < p±/n < 1 GeV/c), all meson and baryon data points fall 
onto an uniform curve. The % data (not shown) are significantly above the V2 of other 
mesons. The large fraction of resonance decay contribution to the % yield is known to 
enhance the V2 of % at a given pj_ [49, 50]. 
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FIGURE 7. Azimuthal angular anisotropy V2/n as a function of p± jn for identified particles where n is 
the number of constituent quarks. The line is a polynomial fit to the data points excluding the % data. The 
bottom panel shows the (multiplicative) deviation from the fit line. 
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For p±/n < 0.6 GeV, there appears a small residual but systematic particle dependent 
deviation from the fit curve, which agrees best with the K s results. The residual mass 
and/or strangeness and/or meson-baryon dependence of this deviation can have many 
causes involving microscopic transport and/or collective flow phenomena. Its detailed 
understanding could play an important role in the demonstration of physics processes 
that lead to the constituent quark level azimuthal flow at hadronization. 

We believe that the particle dependence on vi, which is largley explained by n-scaling, 
requires a V2 distribution at the constituent quark level. The scaled azimuthal angular 
anisotropy (v2/n) may be interpreted as the constituent quark anisotropy just prior to 
the hadron formation. The n-scaling works both for strange and nonstrange quarks, in- 
dicating that the azimuthal flow anisotropy is the same for all three 'light' flavors, and 
by extension, that the collective quark flow is the same for the three flavors. The preci- 
sion of the n scaling leaves very little space for the presence of gluon (fragmentation) 
participation in the intermediate momentum particle production process. This suggests 
that semi-hard gluons have effectively disapperad as independent partonic degrees of 
freedom in the final hadron formation. 

We have shown in this section that the hadron formation dynamics, at RHIC at 
0.6 < p_j_/n < 2, GeV/c is very different from the parton fragmentation picture where the 
leading parton plays a dominant role in determining properties of the final state hadron. 
The measured features for these intermediate hadrons produced at RHIC require that 
all quark ingredients (n = 2 for mesons and n = 3 for baryons) play an approximately 
equal role in the hadron formation and that the hadron properties are determined by the 
sum of contributing partons. 

3. DECONFINED BULK MATTER 

3.1. Bulk properties as function of centrality 

The global study of soft hadron yields and spectra implies that we can fully charac- 
terize the properties of bulk matter at time of hadronization by adding up strangeness, 
entropy, etc., contained in the final state hadrons [51]. This is done particularly eas- 
ily within the statistical hadronization scheme, which have been successfully applied to 
describe stable particle production. 

As a first step, this leads to a quantitative understanding of the variation of the bulk 
properties with centrality [52]. In such an analysis of ^/snn = 200 GeV STAR and 
PHENIX results, the contents in entropy, baryon number, strangeness, and energy grows 
linearly with the participant number for A > 20. This implies that the conversion of 
colliding nuclear matter into bulk parton matter is at RHIC rather independent of the 
size of nuclear overlap and complete, including all matter. In our above discussion of V2, 
we have tacitly assumed that the properties of bulk partonic matter formed in non central 
collisions are not significantly dependent on centrality, as the present study confirms 

In most central 5% of the reactions, one finds that the net baryon density per unit of 
rapidity is at mid-rapidity d(B — B)/dy = 14±2, the strangeness yield ds/dy = 135 ± 10 
and the entropy yield dS/dy = 4900 ±400. The errors comprise the uncertainty of the 
data and the chemical scheme used. The strangeness per entropy yield s/S = (2.9 ± 
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0.3) 10 3 is more than 4 times enhanced compared to the AGS energy scale. For the 
most peripheral reactions, this ratio is s/S — (1.9 ±0.3) 10 -3 which shows the influence 
of the fireball expansion dynamics on production of strangeness. 

Given the experimental yields of particles at central rapidity, the intensive properties 
of bulk matter at hadronization, such as P pressure, £ energy density, o entropy density, 
E/TS = £/T<7, seen in figure 8, follow summing the properties of individual particle 
fractions. They are given for the two chemical models in figure 8 on left. The lines 
guide the eye, the actual results are the squares and circles centered at the mean trigger 
centrality of the data considered. These results are for the full chemical non-equilibrium 
(squares) and for strangeness chemical non-equilibrium (circles) respectively. 
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FIGURE 8. Left: pressure P, energy density e = E/V, entropy density S/V and E/TS and; Right: 
temperature T, light quark phase space occupancy y q , the ratio of strange to light quark phase space 
occupancies y s j y q and the chemical potentials (B for baryochemical jUg and S for strangeness fig) as a 
function of centrality. After Ref. [52]. 

To obtain the physical properties, a fit of the statistical hadronization model parame- 
ters had been performed and the results are shown on right: the hadronization temper- 
ature T, phase space occupancies Y q ,Y s /Yq, and the chemical potential /1b,s of baryon 
number and strangeness, respectively. We see the expected rise of strangeness occu- 
pancy yield with size (centrality) of the system, related to the extended lifespan dur- 
ing which strangeness is produced [33, 34]. Otherwise there is remarkable stability of 
the bulk properties with centrality of the collision. The bulk matter created at RHIC 
at v/^nn = 200 GeV, small or large in size and independently of models of chemical 
(non)equilibrium has a baryochemical potential /Ib = 25 ± 1 MeV. 
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3.2. The Phase Boundary 



The hadronization temperature T = 140-155 ± 8 MeV is near but below the phase 
transformation boundary for QCD with 2+1 flavors. The infinite matter static equilib- 
rium transformation is expected at T = 164 ±4 MeV, for more detail about the structure 
of the cross-over region, see figure 9. 
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FIGURE 9. The phase diagram at low baryon density in physical units, results of Ref. [53]. Phase 
crossover ends in the square which shows the endpoint where 1st order phase transition sets in. 

The lower RHIC freeze-out temperature we reported above is associated with a rapid 
expansion of the RHIC fireball: the flow of color charge pushes the vacuum and the bulk 
breakup occurs suddenly from an over expanded supercooled condition. Quantitatively 
the supercooling is of the required magnitude [47]. The sudden hadronization mecha- 
nism also explains why single freeze-out analysis of particle spectra at SPS and RHIC 
succeed [54, 55]. 



4. HEAVY QUARKS 
4.1. Theoretical remarks 

Heavy quarks (charm, in future bottom) are produced mostly in initial parton scatter- 
ing (gluon fusion) and in any case during the initial primary and/or very high temper- 
ature phase of the collision. Therefore, heavy quark measurement can probe the initial 
parton flux, the dynamical evolution and the quark energy loss in dense medium. If 
heavy quarks are found to participate in the collective motion of the medium (radial or 
elliptic flow), this will lend further confirmation for parton collectivity. 

At RHIC, one expects that multiple heavy quark-antiquark pairs will be produced 
in an individual AA collision. Therefore one can expect a novel mechanism of 
formation by recombination [56]. Then, there should be a contribution to heavy quarko- 
nium formation which utilizes combinations of initially uncorrected quark and anti- 
quark, leading to a quadratic increase with the total number of heavy flavor quarks in 
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the event. This can occur within the deconfined phase, considering that recent Lattice 
QCD (LQCD) calculations indicate that spectral functions of pseudo scalar and vector 
mesons have non-trivial shapes at a temperature above the critical T c [57]. In particular, 
heavy quarkonium such as J/ty/ may survive at a temperature above l.6T c [58, 59]. 

Charm quark transport dynamics in dense nuclear medium will provide unique probes 
to the QCD properties of the medium. If the initial temperature is very high, T ~ 500 
MeV or higher, the yield of total charm quarks can also be increased through thermal 
gluon-gluon scatterings. Possible suppression of charm mesons at high p^ will test 
the energy loss dynamics of charm quark propagation in a QCD medium. Theoretical 
calculations have predicted a reduced medium induced energy loss for heavy quarks 
and the high p± suppression of charmed hadrons should not be as strong as light 
hadrons [60, 61]. For a recent review of theoretical charm quark situation see Ref. [62]. 



4.2. RHIC experimental status 

Both the STAR and PHENIX collaborations at RHIC have been pursuing vigorous 
heavy quark physics programs, both in analysis of current data and in future detector 
upgrade plans, to provide better capabilities for heavy quark measurements. 

The total charm quark pair production cross section (a C c) is an important constraint on 
the collision dynamics and the heavy quark evolution. Both STAR and PHENIX have 
presented results on the o C c measurement of pp collisions from charm semi-leptonic 
decays. In addition, STAR has also derived an equivalent o C c for pp collisions based 
on direct reconstruction of hadronic decays of charm mesons from dAu collisions. 
Figure 10 shows the p^-spectra of electron and D° from STAR. The PHENIX pre- 
liminary non-photonic electron data are represented by the fitted line with a reported 
measurement o C c = 709 ± 85(stat)+332 — 281(sys) \lb [63]. STAR has measured o C c — 
1300 ± 200 ± 400 lib and a mean transverse momentum for D° of (p ± ) = 1.32 ± 0.08 
GeV/c from direct D° reconstruction [64] . A next-to-leading order pQCD calculation of 
the charm quark production cross section [65] has yielded o C c = 300 — 450 lib, sig- 
nificantly below the STAR measurement and at the lower end of the PHENIX range of 
uncertainty. 

Several comments on the cross section measurements are in order. When measuring 
the charm cross section through semi-leptonic decay, the quality of the electron data for 
p_L <, 1 GeV/c is significantly deteriorated because of a large combinatorial background. 
An electron of p^ ~ 1 GeV/c typically comes from the decay of a D meson with p^ ~ 2 
GeV/c, which is significantly beyond the average p± of 1.32 ±0.08 GeV/c reported by 
STAR. Therefore, one has to extrapolate to the low pj_ region by over a factor of two 
to obtain the total charm cross section. Such an extrapolation is often model dependent 
and has a large uncertainty. The semi-leptonic decay branching ratios for D°, D* , D ± 
and Ds are different. The electron yield from decays of these D mesons depends on the 
relative yield which is another important contribution to the uncertainties of the charm 
production cross section derived from electron measurement. The direct reconstruction 
of the D decay kinematics provides a broad coverage of p^ and does not suffer from 
the same uncertainties as the leptonic decay electron measurement. However, present 
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FIGURE 10. /?x distribution of D° mesons and non-photonic electrons from semi-leptonic decays of 
charm mesons. 



STAR measurements of D meson yields using event-mixing methods from TPC tracks 
suffer from limited statistics. A future vertex detector upgrade capable of measuring 
the D decay vertex displacement is essential for both STAR and PHENIX heavy flavor 
physics programs. 

Figure 11 presents the STAR preliminary transverse momentum spectrum of D 
mesons from dAu collisions normalized by the number of binary collisions, where the 
p± shapes of D* and D ± are assumed to be the same as that of D° [66]. The shape of 
the p± distribution coincides with the bare charm quark p± distribution from the Fixed- 
Order-Next-Leading-Log (FONLL) pQCD calculation from M. Cacciari et al. [67]. If 
a fragmentation function such as the Lund fragmentation scheme [68] or the Peterson 
function [69] is introduced for D meson production, the resulting p± distribution will 
be significantly below the measurement at the high p± region. This observation raises 
an outstanding question regarding the p± distribution and the formation mechanism of 
D mesons in hadro-production. Recently, a k± factorization scheme has been found to 
significantly change the D meson pj distribution from nuclear collisions as well [70]. 

The fact that the D meson p± distribution can be better described by the p± of bare 
charm quarks from pQCD calculations has been observed in previous fixed target exper- 
iments [71]. With a fragmentation function such as the Peterson function the D meson 
p± distribution is too steep to explain the measured p^ distribution. In order to match 
the calculation with the data one has to introduce a k± kick to the parton distribution. 
The scale of (kjj is on the order of 1 (GeV/c) 2 , much larger than the typical Aqcd scale 
for strong interaction. Furthermore, the Feynman xf of the D meson distribution was 
also found to coincide with the bare charm quark xp distribution [71]. The fragmenta- 
tion function will have a large impact on the xp distribution from bare charm quark to D 
meson, which cannot be negated by introducing any kp longitudinal boost of reasonable 
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FIGURE 11. STAR preliminary measurement of p± distribution of D mesons from dAu collisions 
normalized by the number of binary collisions. The shape of the p± distribution is compared to pQCD 
FONLL calculations. 



scale as in the case for kick in the transverse momentum direction. 

The transverse momentum distribution of particles produced at RHIC energies are 
considerably flatter than those at lower energies. The k± kick scheme does not change 
the shape of the p± distribution significantly. The STAR measurement of the D meson 
p± distribution suggests either that the charm quark fragmentation may be close to a 
delta function or that a charm formation mechanism such as the recombination model 
may be important in hadro-production. The recombination model takes a charm quark 
and combines it with a light quark, presumably of low p±, to form a charmed meson. 
Therefore, the p± of the meson is not significantly different from that of the bare charm 
quark. In that way, measurements of charmed meson production provide unique probes 
for the hadron formation dynamics and for the transport dynamics of heavy quarks in 
the dense nuclear medium produced at RHIC energies. 

Observation of heavy quark hydrodynamic flow would indicate that heavy quarks, 
once created in the initial state, must have participated in the partonic hydrodynamic 
evolution over a sufficiently long period of time to reach a substantial flow magnitude. 
This would be a unique probe for the early stage of a partonic phase [72]. Figure 12 
shows preliminary STAR [73] and PHENIX [74] measurements of the V2 for electrons 
from charm leptonic decays, which have been demonstrated to be closely correlated with 
charmed meson V2 [50] . 



5. SUMMARY 



Among the most important aspect of many results presented is that the gluon degrees of 
freedom are not explicitly manifested in the hadron formation. Quarks are the dominant 



Hadronization and Quark Probes of Deconfinement at RHIC 



February 2, 2008 



15 



0.3 
0.25 

0.2 
0.15 

CM 
> 0.1 

0.05 



-0.05 
-0.1 



■^■STAR preliminary 
0-80% central 



Tion-photonic (e + +e")/2 



• PHENIX prelim. 
Min.Bias 



0.5 



1.5 



v 2c-quark _ v 2q 
v 2c-quark = " 



I 



p T (GeV/c) 



2.5 




FIGURE 12. Preliminary measurement of non-photonic electron V2 as a function of pj from PHENIX 
and STAR. The curves are from a calculation by V. Greco et al. [72]. 



degrees of freedom at the boundary of quark-hadron transition. It is a critical step to 
firmly connect this mainly experimental insights on the properties of the quark matter 
at the boundary of hadronization with the LQCD calculations. The disappearance of the 
gluon degrees of freedom from the initial state and the emergence of constituent quarks 
at hadronization are some of the critical conceptual questions to be addressed in the 
imminent future [75]. 

Strange (anti) baryon production enhancement at RHIC as function of centrality, 
energy and p^ is significant. For <, 5 GeV/c, relative baryon to meson yields have 
shown distinct features that are drastically different from the fragmentation processes 
in elementary collisions. A novel feature of meson and baryon dependence has been 
observed in the nuclear modification factor and the angular anisotropy V2 of n, K ± , 
Ks, K* , p, A and S particles at the intermediate p± of 2-5 GeV/c. A constituent 
quark number scaling has been observed for the V2 measurement. These experimental 
measurements suggest that quarks have developed a collective V2 as a function of p±; 
and the hadron formation at the intermediate p± is likely through multi-parton dynamics 
such as recombination or coalescence process. 

This physical picture emerging from the experimental measurements complements 
the Lattice QCD results. Spectral function calculations have indicated that hadrons, 
particularly heavy quarkonium, do not melt completely at critical temperature. It appears 
plausible that the constituent quark degrees of freedom or even hadron-like quasi- 
particles play a dominant role at the hadronization of bulk partonic matter though further 
confirmation of the picture from LQCD is needed. We note the related theoretical 
models [76], which invoke the notion of quasi-hadrons to describe properties of the 
dense matter. 
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We have also described how heavy (charm) quark production and its transport dy- 
namics in dense nuclear medium probe QCD properties of the matter. The charm quark 
flow measurement will provide a significant insight on recombination or coalescence 
hadronization mechanism and partonic collectivity of the dense matter. Charmonium 
enhancement arises from a dense charm rich state. Future detector upgrades from STAR 
and PHENIX will greatly enhance their heavy quark measurement capabilities at RHIC. 

Despite intriguing experimental observations of hadronization from a deconfined 
bulk partonic matter, a smoking colt signatures for the quark-hadron phase transition 
remain elusive. On the other hand, the totality of experimental results can be understood 
invoking deconfinement. 
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